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I.  INTRODUCTION 


Spall  fragments  produced  behind  armor  during  impact  must  be 
characterized  in  detail  so  that  damage  potentials  of  shaped 
charges  and  other  penetrating  devices  can  be  assessed 
adequately.  The  characterization  usually  consists  of  mass 
distributions,  directions  of  travel  and  individual  fragment 
velocity  or  speed.  Such  descriptions  have  traditionally  been 
obtained  by  collecting  fragments  after  each  field  test  through 
the  use  of  a  relatively  soft  recovery  medium  such  as  wall  board. 
At  present,  another  technique,  time  sequenced  flash  x-rays  of 
spall  just  after  formation  is  being  developed  to  complement  or 
replace  the  recovery  (collection)  method,  but  until  improved 
further,  the  recovery  technique  is  apt  to  remain  in  use.  The 
recovery  method  is  not  without  inherent  limitations, 
particularly  with  respect  to  the  ability  to  accurately  determine 
the  initial  spall  fragment  velocity  before  impact  with  the 
recovery  material.  Thus,  it  is  necessary  to  examine  the 
technique  from  time  to  time  for  possible  improvement. 

The  velocity  determination  is  obtained  through  use  of  an 
empirical  equation  calibrated  for  the  recovery  material 
utilized.  Generally  such  equations  provide  a  relationship 
between  the  striking  velocity  and  penetration  depth  for  given 
fragments  with  known  masses.  The  presented  areas  of  the 
fragments  are  a  variable  of  the  system  and  in  practice  are 
either  measured  or  assumed.  Simulated  fragments  such  as 
spheres,  cubes  or  short  cylinders  are  often  used  and  in  these 
cases  the  data  obtained  is  rather  uniform,  but  does  include  some 
scatter  due  to  uncertainty  in  fragment  orientation  and 
variations  in  recovery  material  properties.  The  application  of 
these  equations  to  irregular  spall  fragments  gives  rise  to 
considerably  more  uncertainty  since  the  individual  presented 
areas  differ  significantly.  Further,  presented  areas  associated 
with  the  extreme  orientations  create  an  even  wider  dispersion 
which  is  not  easily  taken  into  account. 

Zook*  includes  a  historical  summary  of  such  calibration 
equations.  The  most  widely  used  calibration  expressions  are 

,  2 

perhaps  those  determined  empirically  by  Project  Thor  , 


1.  Zook,  John.  "An  analytical  model  of  Kinetic  Energy 
Projectile/Fragment  Penetration."  BRLMR  2797,  October  1977 

2.  Ballistic  Analysis  Laboratory,  Johns  Hopkins  University. 
Project  Thor  Technical  Report  No.  50:  "The  Calibration  of  a 
Collection  Medium  for  the  Determination  of  Particle  Velocity." 
July  1962. 
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Whiteford  and  Regan  ,  and  Collins  .  The  mathematical 

a  h  c 

expression,  V  *  KM  P  A  ,  appears  in  their  work  together  with  an 

evaulation  of  the  constants,  K,  a,  b,  and  c  fitted  to  various 
sets  of  data  which  reflect  a  number  of  wall  board  materials  and 
fragment  characteristics. 

In  the  expression,  V  is  the  striking  velocity  and  M,  P  and  A 
are  respectively  the  mass,  penetration  depth,  and  the  average 
cross  sectional  area  of  the  fragment  type  utilized.  The 
inclusion  of  the  variable  A  would  presumably  adjust  the 
equation  adequately  for  various  fragment  shapes  (various 
average  presented  areas).  However,  an  examination  of  various 
sets  of  constants  show  that  they  differ  significantly  as  does 
the  predicted  velocity  for  a  given  set  of  fragment/penetration 
c ond i t i on s  . 


II.  THE  MODEL 

This  study  begins  with  the  following  hypothesis;  The 
resistance  force  acting  on  a  projectile  during  penetration  of  a 
thick  impact  medium  (thick  enough  to  bring  the  projectile  to 
rest)  is  not  constant.  Near  the  beginning  of  penetration, 
where  the  velocity  is  high,  the  force  is  high  and  rapidly 
changing,  similiar  to  a  hydrodynamic  drag  force 


At  modest  velocities,  the  force  is  derived  from  viscous 
effects : 


Near  the  end  of  penetration,  where  the  velocity  is  low, 
the  force  is  primarily  frictional,  and  essentially  constant: 


Apparently  the  exponent  of  v  varies  over  the  course  of  the 
penetration.  In  terms  of  the  dependent  variable,  one  can  say 
the  exponent  of  v  varies  over  the  range  of  v.  Therefore,  the 
exponent  of  v  is  itself  a  function  of  v,  varying  from  nearly  2 
at  the  beginning  of  penetration  to  nearly  0  at  the  end. 


.  Whiteford,  C.  W.  and  Regan,  J.  M.,  "The  Determination  of 
the  Striking  Velocity  of  Steel  Fragments  by  their  Mass  and 
Penetration  into  Witness  Material",  BRLMR-1333,  April  1961. 

4.  Collins,  John  A.  "Fiberboard  Calibration  for  Determination 
of  Fragment  Velocities",  Eglin  Air  Force  Base  Technical  Report 
AFATL-TR-73-193 ,  September  1973. 
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Algebraically,  then,  the  assumed  basis  of  this  model  is: 
f  =  -kv  g(v) 

In  this  force  law  g  is  some  function  of  v  for  which 
o<g<2  for  0<v<vo 

where  vq  is  the  striking  velocity.  This  is  all  that  is  known 

about  the  function  g.  Since  most  functions  can  be  represented 
by  a  polynomial  (e.g.  Taylor  series),  this  is  a  logical  choice 
for  g.  Specifically,  a  single  term  "polynomial"  is  chosen  for 
simplicity.  Thus  the  assumed  velocity  dependence  is: 

g ( v  )  =  ava 

In  the  force  law  k  =  k(v)  to  the  extent  of  reconciling 
dimensions  (i.e.  numerically  constant,  dimensionally  variable). 
Following  conventional  theory,  the  constant  depends  on  the 
density  of  the  medium  being  penetrated  (p)  and  the  presented 
area  of  the  projectile  (A).  Finally, 

av  a 

f  =  - c  p  Av 

where  c  is  a  calibration  constant  for  the  impact  medium  (i.e. 
determined  by  calibration  experiments),  P  is  the  density  of  the 
Impact  medium,  A  is  the  presented  area  of  the  projectile,  and  o 
and  a  are  as  yet  unknown  constants.  They  must,  however,  be 
positive  to  fit  the  model  presented.  Note  that  at  this  point 
conventional  dimensional  analysis  will  yield  the  familiar 

2 

hydrodynamic  drag  force,  CpP Av  . 

Since  the  exponent  must  be  dimensionless,  let 

*  =  v/va 

O 

where  V  is  a  dimensionless  constant  (as  yet  unknown).  Thus: 


Inspection  of  Equation  1  reveals  that  *.' =  2  and  a  i  0  is 
required  if  the  exponent  is  to  vary  from  0  to  2.  The  polynomial 
form  of  the  approximation  of  the  exponent  is  not  valid  at  v  =  0, 
since  it  leads  to  a  non-zero  force  there.  In  lieu  of  choosing 

g 

a  two  term  polynomial  for  the  exponent,  (e.g.  av  +  (3  )  the  range 
of  v  will  be  limited  to  v/v0  >  0.  This  will  cause  no  inconven¬ 
ience  in  the  analysis  of  actual  experimental  data,  nor  will  it 
result  in  any  appreciable  error. 
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From  Newton's  second  law,  if  the  recovery  medium  causes  no 
erosion  of  the  projectile: 


M 


dv 

at 


c 


.  2 (v/v  )' 
•  p cAv  o 


(2) 


where  M  is  the  mass  of  the  projectile.  Now  the  derivative  may 
be  rewritten  as  follows.  Since  v  =  v  (x,t). 


dv  =  dx  ♦  1^-  dt 
ax  ct 


dv  _  0V.  dx  +  22. 
dt  ~  ax  dt  at 

but 


dv 

at 


0 


dv  _  dv  dx 
dt  dx  dt 

where  dx/dt  is  simply  v.  Substituting  into  Equation  2  will 
yield  the  velocity  as  a  function  of  penetration  depth,  x. 


Mv  *  -fcAv2Cv/vo)3 
Separating: 


vl-2(v/v0)' 


dv  -  “ 


pcA 


dx 


(3) 


Both  sides  are  integrated  over  the  course  of  the  penetration: 
JCA  t'  dx  .  f‘  dv 

"  M  Jo  Jv 

0 

where  P  is  the  penetration  of  the  projectile  at  v  -+  0, 
approximately  equal  to  the  total  penetration. 

v 

*>cA  _  (  0  1-2 (v/v  )a  . 

P  ■  ^  v  v  oJ  dv 


III.  APPLICATION  OF  THE  MODEL:  METHOD  I 

Equation  contains  two  calibration  constants,  a  and  c.  To 
evaluate  these  two  unknown  empirical  constants,  recall  that  c  is 
a  constant  associated  with  the  impact  medium,  and  is  not  a 

I 
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function  of  projectile  area,  mass  or  velocity  (assumed ly ) .  Thus 
the  second  "equation"  needed  to  evaluate  the  two  unknowns  a  and 
c  is  c  *  constan'.,  and  the  value  of  a  is  sought  which  best 
satisfies  thi"  condition. 

Collins'  impact  medium  calibration  using  57  mg  steel 
spheres  (Table  A-l)  is  used  as  an  example.  Equation  4  is  solved 
for  c  . 


c 


M 

PAP 


Kv0) 


(5) 


where 

I(v0)  ■  jV°  v'-!WV‘  dv  <6> 

represents  the  integral  on  the  right  in  equation  4.  In  Collins' 
test,  the  impact  medium  was  Celotex,  with  a  density  of 


P=  272  kg/m3  (7) 

For  57  mg  spheres  , 

M  -  5. 7  X  10-5  kg  (8) 

The  presented  area  at  impact  is  constant  for  all  orien¬ 
tations,  and  is  found  from  the  sphere's  mass: 

A  -  4.54  X  10'6  m2  (9) 

R  *  1.20  X  10_3m  (10) 


Substituting  equations  7,  8,  and  9  into  Equation  5: 


Uv  ) 

c  =  .04616  -=-2- 


(11) 


Collins'  calibration  data  (Table  A-l)  gives  v  and  P  for 

o 

numerous  shots  of  the  steel  spheres  into  the  Celotex  impact 
medium.  Tests  yielding  penetrations  of  less  than  .051m  showed 
erratic  results,  and  so  were  disregarded  in  this  study.  The 
following  method  was  used  to  evaluate  a  and  c:  a  is  a 
parameter,  and  is  assigned  various  values.  For  a  given  data 
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point  (P,  v^),  I  is  found  by  numerical  integration  of  Equation  6. 

The  trapezoidal  rule  with  1000  uniform  intervals  was  used  for 
accuracy  consistent  with  the  accuracy  of  the  data.  The  value 
for  I  and  the  corresponding  value  for  P  (for  the  given  data 
point)  is  then  substituted  into  equation  11  to  determine  c. 

This  procedure  is  repeated  for  each  data  point  (P,  v  )  and 

the  mean  Cc)  and  standard  deviation  (<r)  of  the  resulting  values 
of  c  are  computed.  A  new  value  for  the  parameter  a  is  then 
assumed,  and  the  whole  process  is  repeated.  Since  c  should  be 
constant,  a  should  be  zero.  Of  course,  this  is  not  reasonable 
in  an  experiment.  Instead,  a  minimum  is  sought.  c  is 
normalized  by  the  mean  (c")  and  a/c  is  plotted  versus  the 
parameter  a.  This  is  done  in  Figure  1.  A  clear  minimum  exists. 
Thus  we  conclude  that  the  proper  value  of  the  parameter  a  is 

a  *=  0.20  (12) 

For  this  choice  of  a,  Collins'  57mg  sphere  calibration  data 
leads  to  an  average  value  for  c  of 

c  »  50.3  (13) 

The  computations  are  summarized  in  Table  1,  and  were 
obtained  from  Equation  5  with  numerical  integration  of  Equation 
6.  We  are  now  ready  to  use  Equation  4  to  predict  the  velocity 
of  fragments  impacting  the  Celotex  impact  medium.  Collins' 
data  for  small  fragments  (approximately  the  same  mass  as  the 
calibration  spheres),  presented  in  Table  A-2,  will  be  studied. 

First  Equation  5  is  rewritten 


I(vo) 


(  14) 


Now  I (v  )  is  defined  by  equations  6  and  12.  I(v  )  can  be 
•  o  0 

numerically  integrated  and  plotted  as  a  function  of  vq.  This  is 

done  in  Table  2  and  Figure  2;  the  relationship  turns  out  to  be 
nearly  linear.  This  simplifies  the  prediction  of  fragment 
velocities,  since  a  graphical  solution  is  not  required.  For 
velocities  greater  than  lOOm/s,  the  equation  of  the  curve  in 
Figure  2  is: 


I(v0) 


(15) 
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FH'llRi:  1.  IViramct ric  Variation  of  Calibration  Data. 


Table 


1.  Computation  of  Impact  Medium  Calibration  Constant  c. 


1606. 

.155 

1602. 

.155 

1600. 

.155 

1588. 

.148 

1567  . 

.151 

1535  . 

.138 

1529. 

.151 

1479. 

.135 

1375  . 

.133 

1369. 

.138 

1362. 

.135 

1354. 

.135 

1351 . 

.  141 

1340. 

.135 

1331 . 

.135 

1023  . 

.106 

1015. 

.106 

1007  . 

.099 

995  . 

.104 

993  . 

.101 

988. 

.101 

978  . 

.106 

969  . 

.098 

962  . 

.098 

962  . 

.098 

787  . 

.077 

753  . 

.077 

693. 

.077 

691  . 

.071 

672  . 

.064 

663  . 

.071 

652  . 

.071 

652  . 

.064 

621  . 

.064 

617. 

.064 

615. 

.057 

612  . 

.064 

565  . 

.057 

c 

I C  v  ) 

0 

50.10 

168.23 

49.98 

167.83 

49.92 

167.64 

51.92 

166.46 

50.26 

164.40 

53.94 

161.27 

49.12 

160.68 

53.27 

155.79 

50.53 

145.60 

48.51 

145.02 

49.35 

144.33 

49.08 

143.55 

46.90 

143.25 

48.61 

142.18 

48.31 

141.30 

48.39 

111.13 

48.05 

110.35 

51.09 

109.56 

48.11 

108.39 

49  .45 

108.19 

49.22 

107.70 

46.47 

106.72 

49.85 

105.83 

49.53 

105.15 

49.53 

105.15 

52.71 

87.93 

50.70 

84.58 

47.14 

78.64 

51.00 

78.44 

55.22 

76.56 

49.19 

75.67 

48.48 

74.58 

53.79 

74.58 

51.56 

71.49 

51.28 

71.10 

57.41 

70.90 

50.92 

70.60 

53.37 

65.91 
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Table  2.  Numerical  integration  of  the  integral 


vo 

INTEGRAL 

100. 

16.90 

200  . 

28.27 

300. 

38.93 

400. 

49.26 

500. 

59.39 

600. 

69.40 

700  . 

79.34 

800  . 

89.22 

900  . 

99.06 

1000. 

106.88 

1100  . 

118.68 

1200. 

128.47 

1300  . 

138.36 

1400. 

148.05 

1500  . 

157.84 

1600  . 

167.64 

1700  . 

177.44 

1800. 

187.26 

1900. 

197.08 

2000. 

206.91 

Thus  equation  4  can  be  rewritten 


pcA 

f 


P 
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Despite  the  complexity  of  the  original  model,  the  depth  of 
penetration  of  the  projectile  is  simply  directly  proportional 
to  the  striking  velocity  (this  may  not  be  the  case  for  other 
sets  of  data).  Rearranging, 


IOpcA  p 
o  M 


80 


(16) 


o 

Substituting  p  *  272  kg/m  and  c  ■  50.3 

v  =  1.37x10s  ^  P  -  80 
o  M 


(17) 


This  can  be  written  in  the  general  form 


(18) 


where  P  is  the  slope  and  b  the  ordinate  intercept  of  the  line. 
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Collins’  data  for  small  fragments  (see  Table  A-2)  is 

plotted  in  Figure  3.  Fragments  for  which  v  was  less  than 

o 

lOOm/s  were  disregarded  (total  of  three).  A  linear  regression 
of  this  data  was  performed,  with  the  following  results: 

slope  (/?)  *=  1.14  X  105  kg/m2  s 
ordinate  intercept  (b)  *  143  m/s 

This  line  is  plotted  in  Figure  3  as  the  solid  (actual) 
line.  Also  plotted  is  the  dashed  (predicted)  line,  which  is 
based  on 

slope  (/?)  -  1  .  37  X  1  05  kg/m2  s 
ordinate  intercept  (b)  *  -80  m/s 

from  Equation  17  above.  A  close  look  at  Figure  3  reveals  that 
for  impact  velocities  between  1000  and  1500  m/s,  the  described 
method  predicts  within  +_  5%  of  the  true  value  and  within  +_  10% 
between  800  and  2600  m/s.  The  error  at  200  m/s  is  about  100%. 
The  method  is  more  accurate  at  high  impact  velocities  by  virtue 
of  the  exponent  of  v.  Also,  one  would  expect  more  consistent 
penetration  data  at  high  velocities.  Also,  recall  that  Equation 
17  is  not  valid  at  striking  velocities  less  than  100  m/s  (e.g. 
penetration  at  or  near  zero  must  be  excluded  from 
consideration)  . 


IV.  SIMPLIFYING  THE  MODEL 


Consider  the  velocity  of  the  projectile  as  it  penetrates 

the  recovery  medium.  Equation  3  is  integrated  from  0  to  x  and 

from  v  to  v(x).  a  *  0.2  from  Collins'  data  was  used  for 
o 

definiteness,  although  a  ■  1  was  also  tried  with  similar 
results. 


-pc 


vl-2(v/vo)0,2 


dv 


(19) 


To  simplify  notation,  the  integrand  will  be  referred  to  as 
f(v).  Solving  for  x, 


x 


1_  M 
Pc  A 


f(v)  dv 


(20) 


11 


4. 


FIGURF.  3.  Results  of  Mr  rhod 


Reducing  the  equation  to  the  customary  form  v  =  v(x)  is 

not  possible,  so  the  following  discussion  will  be  somewhat 
awkward.  The  integrand  f(v)  is  plotted  in  Figure  4.  The 
independent  variable  is  of  necessity  v,  the  instantaneous 
projectile  velocity.  For  any  striking  velocity  vq  and 

instantaneous  projectile  velocity  v,  the  area  under  the 

appropriate  v  curve  between  v  and  v  (as  a  fraction  of  the 
o  o 

total  area  under  the  vq  curve)  is  proportional  to  the  depth  in 

the  impact  medium  at  which  the  velocity  v  is  reached  (as  a 
fraction  of  the  total  depth  penetrated). 

Figure  4  thus  gives  important  insight  into  the  model.  The 
projectile  slows  down  extremely  rapidly  in  the  beginning  of  its 
penetration  into  the  medium;  approximately  one-half  of  its 
velocity  is  lost  in  the  first  one-tenth  of  penetration.  Thus 
the  hydrodynamic  regime,  where  the  resistance  force  is 

2 

proportional  to  v  ,  is  over  very  quickly,  with  little 
contribution  to  the  total  penetration.  In  fact,  the  great 
majority  of  the  penetration  occurs  in  a  narrow  band  of  velocity 
exponent,  say  from  1/4  to  3/4.  Thus  assuming  a  constant  value 
for  the  exponent  of  velocity  in  the  resistance  relationship  is 
reasonable.  This  assumption  should  lead  to  acceptable  accuracy 
over  a  sizable  range  of  striking  velocity  provided  the  value  of 
the  exponent  is  chosen  properly.  How  is  this  done? 

To  begin,  the  resistance  force  (Equation  1)  with  a  *  0  will 
result  in  a  constant  for  the  velocity  exponent: 


f  =  -pcAvy 


(21) 


where  y  is  a  constant  to  be  determined.  This  force  is  zero  at 
v  *  0,  so  the  range  of  applicability  of  v  is  0 <  v <  v  .  Equation  4 
now  b e come s : 


To  avoid  the  singularity  at  r  *  2,  Let 


(22) 


y=  2  -  I  (23) 

8  o  <  a  <  oo 

p  turns  out  to  be  a  more  useful  parameter  than  y  . 
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F I  C;UR1:.  4.  Fragment  Penetration  as  a  Function  of  Inst  ancons  Fragment  Velocity 


Sub s  t i tut ing  : 


Integrating  : 


dv 


(  2  A ) 


This  can  be  rewritten 


v 

o 


M-6A6Pe 


(25) 


which  has  the  mathematical  form  extensively  used  for  velocity 
calibration  of  recovery  media  as  indicated  previously.  Note 
that  ft  *  1/2  is  the  solution  for  a  constant  resistive  force, 

0  m  1  represents  a  force  proportional  to  the  first  power  of  the 
velocity,  and  /3  »  1  corresponds  to  the  hydrodynamic  drag  force. 


V.  COMPARISON  TO  EXISTING  CALIBRATION  EXPRESSIONS 

In  this  section  Equation  25  is  compared  to  calibration 
expressions  (having  the  same  form)  as  utilized  by  Project 
Thor (ref  2),  Whiteford  and  Regan(ref  3),  Collins(ref  4),  and 
others.  One  form  frequently  encountered  is 

v  «  K.MaPbAC  (26) 

o  1 

In  cases  where  fragment  area  and  mass  can  be  assumed  to  be 

related  by  A  *  KM2^3  (see  discussion  of  Method  IV  below), 
another  common  expression  results  in 

v  -  K,MgPb  (27) 

o  2 

The  exponents  a,  b,  c,  g,  and  the  coefficients  and  Kj 

are  given  in  calibration  determinations  as  elements  of 
regression  analyses  and  in  general  are  not  constrained  to 
maintain  any  given  relationship  between  themselves.  However, 
the  current  assumed  model  (Equation  25)  indicates  that  a 

definite  relationship  exists  between  exponents,  namely 

- a  ”  b  *  c  (28) 


-g  -  1/3  b 


(29) 
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These  relations  were  tested  by  sets  of  empirically 
determined  exponents  originally  determined  for  the  equations  26 
and  27.  The  work  of  several  investigators  was  used,  and  is 
summarized  in  Table  3.  In  all  cases  the  fragments  considered 
were  steel.  Further  details  of  the  tests  may  be  found  in  the 
cited  references. 

The  comparisons  provided  by  the  table  indicate  that  the 
empirical  ratios  are  close  to  those  predicted  by  the  model.  The 
ratios  appear  to  be  very  close  for  data  representing  spherical, 
cubicle,  or  short cyl indrical  (1/d  =  1)  fragments.  However, 
there  is  a  greater  difference  for  data  obtained  in  tests  where 
fragments  had  more  extreme  shapes  (bomb,  weapon,  and  bomblet 
fragments).  The  departures  are  caused  by  unavoidably  inaccurate 
measurements  of  the  average  presented  area  of  fragments, 
rotation  of  the  fragment  during  penetration  (i.e.  variation  of 
presented  area)  and  other  problems  in  a  real  fragmentation 
experiment.  Affecting  calibration  shots  as  well  as 
fragmentation  experiments  are  such  factors  as  instrumentation 
accuracy  and  collection  medium  homogeneity.  The  fact  that 
Equation  25  (and  hence  Equations  28  and  29)  results  from  an 
approximation  to  the  originally  proposed  model  may  also 
contribute  to  the  differences,  depending  on  the  validity  of  the 
original  model . 

The  favorable  comparisons  obtained  with  the  more  ideal 
shapes  lend  credence  to  the  approach  taken  here.  Hence,  the 
theory  provides  the  basis  for  further  examination  of  the  factors 
involved  in  the  calibration  process. 

Historically,  the  exponents  and  lead  constants  are  applied 
collectively  in  any  given  expression.  Thus  although  they  may 
vary  from  the  theoretical  values  because  of  the  quality  of  fit, 
amount  and  quality  of  data,  etc.,  their  collective  application 
provides  reasonable  predictions  of  velocity  over  the  range  of 
applicability  of  the  calibration  measurements.  However,  if  the 
exponents  are  related  as  prescribed  above,  one  would  expect  a 
simplification  of  the  calibration  process  and  a  better  fit  over 
the  rest  of  the  data. 


VI.  APPLICATIONS  OF  THE  SIMPLIFIED  MODEL 
A.  Calibration 


Equation  25  is  the  starting  point  in  reducing  experimental 
data.  It  is  repeated  here  for  reference. 


v 

o 


M-W 


(25) 
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The  constants  c  and  ft  are  determined  by  calibrating  the 
impact  medium  by  firing  known  fragments  at  known  (measured) 
velocity  into  the  medium.  Spheres  are  used  here  because  of  ease 
of  launching  and  their  known  presented  area  (A  is  independent  of 
the  sphere's  orientation  at  impact). 


To  determine  c  and  /?  from  the  calibration  data,  rewrite 
Equation  25  : 


Now 


substitute 


A 


sphere 


7r  R  and  M 


s  pher  e 


(30) 


(31) 


In  Vq  =  Bln  P  +  b 


(32) 


where 


v 

o 

P 


c 

p 

p 

M 

R 


h  £  P  l\ 

Bln  \4  8  Psr)  (33) 

Striking  Velocity  of  fragment,  m/s 

3 

Mass  density  of  impact  medium  (Celotex)  *  272  kg/m 

3  3 

Mass  density  of  sphere  (steel)  ■  7833  kg/m  (489  lb/ft  ) 
Constant  dependent  on  impact  medium,  (m/s)'/' 

Constant  dependent  on  impact  medium  (dimensionless) 

Depth  of  penetration  of  fragment  into  impact  medium,  m 
Mass  of  fragment,  kg 

Radius  of  the  steel  calibration  sphere,  m 
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P  and  bj  in  Equation  32  are  determined  from  a  linear 
regression  analysis  of  the  calioration  data.  c  is  then 
determined  from: 


c  =  |  y-  R8  exp  (bj/8) 


(34) 


Table  A-l  shows  velocity  and  penetration  of  57  mg  steel 
spheres  into  a  particular  impact  medium.  The  data  is  plotted 
in  Figure  5.  Regression  analysis  indicates: 

P  -  0.672 
bx  -  8.448 

and  in  Equation  34: 


c 


(1.20x10'3)(.672)  exp 


8.448 

.672 


c  -  8924 

These  values  of  P  and  c  are  assumed  to  hold  for  any  shape 
fragment  that  impacts  this  medium. 

Tables  A-3  through  A-5  show  similiar  data  for  other  size 
spheres  (700mg  to  1059mg).  Regression  analysis  of  these  groups 
shows  P  and  b^  to  vary  slightly  over  this  range  of  "fragment" 

mass,  indicating  that  calibration  spheres  should  be  of 
approximately  the  same  size  as  the  fragment  being  studied.  A 
given  test  therefore  may  require  several  sizes  of  calibration 
spheres.  The  model  (Equation  25)  does  not  predict  any  variation 
of  P  and  b^  (Equation  32)  with  mass.  The  existence  of  the 

variation  is  most  likely  due  to  variations  in  measurement  of  one 
or  more  of  the  variables  (M,A,P,or  v).  Table  4  shows  the  extent 
of  the  variations  for  Collins'  data. 

Although  Equation  32  predicts  a  linear  relationship  between 
In  P  and  In  v  ,  when  the  actual  data  is  plotted  (Figures  5,  6, 
and  7),  the  data  is  seen  to  deviate  from  a  straight  line, 
especially  at  high  velocities.  Dividing  the  graphs  into  two 
more  nearly  linear  regions  does  not  make  much  difference  in  the 
"constancy"  of  P  and  b^  over  this  range  of  mass  (see  Table  4). 

Note  also  that  grouping  together  all  sphere  data  is 
apparently  misleading,  as  it  results  in  a  lower  slope  (P)  than 
any  individual  group,  while  not  appreciably  affecting  the 
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Table  4.  Variation  with  Fragment  Mass  of  Calibration  Constants 
ft  and  b  ^ . 


MASS 

mg 

N 

Combined 

0<v<8 

f*  ,  bl 

Separate 

0<v<6 

e  |  h 

Separate 
!  6<v<8 

»  i  bi  i 

! 

57 

55 

.672 

8.45 

i 

.477  1 

i 

i 

7.61 

i  .994 

\ 

» 

9.20 

700 

69 

'  .671 

;  7.70 

.4  74  : 

1 

7.02 

i 

1  .940 

8.01 

1025  | 

1 

101 

.736 

,  7.81 

| 

1 

.333  j 

j 

6.48 

1.070 

8.17 

ALL 

225 

.565 

7.6  9 

i 

i 

.263 

_ L 

6.37 

1 

ordinate  intercept,  b^  (high  velocity  behavior  is  similiar). 

Physically,  spheres  of  different  mass  impacting  the  medium  at 
the  same  velocity  will  bring  different  amounts  of  energy  to  the 
collision,  so  one  would  expect  them  to  penetrate  differently; 
the  more  massive  sphere  penetrating  more  deeply.  Presented 
area,  being  greater  for  the  more  massive  sphere,  will  mitigate 
this  effect.  Thus,  grouping  together  all  sphere  data  is  also 
physically  incorrect. 

B.  Method  II 


The  first  simplification  of  the  basic  model  assumes  only  a 
constant  value  for  the  exponent  of  v.  The  test  must  measure 
mass  ,  presented  area,  and  penetration  into  the  calibrated  medium 
for  each  fragment.  Calibration  of  the  impact  medium  (see  above) 
is  also  required.  This  type  of  data  is  presented  in  Table  A-2, 
which  was  selected  from  Collins,  Appendix  J  ("Small  Fragment 
Data"),  for  fragments  of  approximately  the  same  mass  as  the 
calibration  spheres  (57  mg.)  Collins  measured  the  velocity  of 
these  fragments  electronically,  and  this  data  is  also  presented. 

To  predict  fragment  velocity,  begin  with  Equation  25. 

Rewr i ting: 
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LN  V  (  m/s ) 
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FIGURE  7.  Sphere  Calibration  of  Impact  Medium  Using  1025mg  Spheres. 
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(36) 


where 


/pc\  (37) 

b  =  6  m^-J 

/?  and  c  were  found  from  the  calibration  of  the  impact  medium, 
and  A,  P  and  M  were  measured,  so  v  can  be  computed.  From  the 
sphere  calibration  data  above, 

J 8.  -  0.672 
c  «  8924 

The  calibrated  value  for  b,  from  Equation  37  is: 
b  -  10.15 

The  true  (from  measured  velocity  data)  values  for  /?  and  b 
in  Equation  36  are  found  by  regression  analysis  of  the  velocity 
data  in  Table  A-2: 

/3  -0.755 

b  -  10.60 

The  measured  velocity  data  is  plotted  as  the  dots  in  Figure  8.  The 
solid  line  is  the  result  of  a  linear  regression  analysis  of  this  data:  the 
slope  (/J)  is  .755  and  the  ordinate  intercept  (b)  is  10.60.  Two  other  lines 
are  plotted.  The  dashed  line  marked  "Method  II"  is  the  regression  of  the 
predictions  from  method  II,  using  /9  -  .672  and  b  -  10.15.  The  dashed  line 
marked  "Collins"  is  the  regression  of  velocity  predictions  based  on  the 

3  b  c 

relationship  v  -  K  M  P  A  .  Here  a,  b,  and  c  are  Collins'  values  as  given  in 
Table  3  ("sph  &  bomb").  Also,  K  -  7102  was  obtained  from  the  57  mg  sphere 
calibration  data  (Table  A-l).  For  this  line,  p  -  .879  and  b  -  11.26.  The 
figure  shows  this  method  to  be  comparable  in  accuracy  to  Method  II  at  the 
low  and  middle  velocities,  but  at  high  velocity  Method  II  is  superior. 

C.  Method  III 

Several  additional  approximate  methods  exist.  These 
methods  significantly  reduce  the  time  and  cost  involved  in 
conducting  the  test,  but  of  course  there  is  an  additional 
sacrifice  in  accuracy. 

The  first  of  these  methods  avoids  the  need  to  measure  mass 
and  presented  area  for  each  fragment  by  introduction  of  a  shape 
factor,  defined  by  a  statistical  sampling  of  the  fragment 
popul a  t ion . 
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(METRIC  UNITS) 


Define  an  average  shape  factor,  k;  assumed  valid  for  the 
entire  fragment  population: 

A  *  kMr  (38) 

k  may  be  found  by  measuring  a  statistically  valid  sample  of  the 
entire  fragment  population.  If  the  sample  is  indeed  valid  this 
method  can  by  very  nearly  as  accurate  as  method  II,  as  the 
example  below  will  show.  Note  the  subtle  difference  in 
the  two  methods:  using  A^  and  for  each  individual  data 

point  in  method  II  was  tantamount  to  using  individual  shape 
factors.  Method  III,  of  course,  uses  one  average  shape  factor 
for  all  data  points.  From  Equation  38: 

In  A  *  r  In  M  +  In  k  (39) 

Thus  a  linear  regression  analysis  of  fragment  presented  area  and 
mass  data  will  yield  values  for  the  constants  r  and  k.  The 
value  of  r  thus  determined  is  assumed  valid  for  each  individual 
fragment;  it  is  a  function  of  the  entire  population  (class)  of 
fragments . 

Mass  and  area  data  for  the  small  fragments  under  study  are 
included  in  Table  A-6.  Using  equation  39,  a  regression  analysis 
of  the  data  yields: 

r  *  0 . 44  , 

k  *  5.56  x  10  (In  k  *  -7.50) 

This  was  repeated  for  all  fragment  data  (Table  A— 6  and  Figure  9) 

resulting  in  r  *=  0.42  and  k  *  4.69  x  10  .  In  this  case  the 
"statistically  valid  sample"  is  the  entire  population.  In  fact, 
the  entire  population  breaks  down  as  follows: 


Ma  8  s 

(mg ) 

n 

r 

k 

-4 

36 

-  80 

51 

.44 

5.56  x  10 

.  „ ~4 

81 

-  104 

50 

.30 

1.44  x  10 

.  »-4 

105 

-  205 

50 

.40 

3.79  x  10 

.  -4 

ALL 

151 

.42 

4.69  x  10 

Relationship  between  Fragment  Area  and  Fragment  Mass  for  the  F.ntire 
Fragment  Population. 


A  similar  analysis  can  be  done  to  determine  the  degree  of 
dependence  of  r  and  k  on  fragment  area. 

To  utilize  this  method  III,  substitute  Equation  38  into 
Equa t ion  2  5  : 


(41) 


Where  b  is  defined  as  Equation  37.  From  the  sphere  calibration 
data  we  get  0  and  c,  while  r  and  k  were  found  from  the  actual 
fragment  data.  So  according  to  the  method  presented  here, 
and  b  are  unchanged: 

0  -  0.672 

b  -  10.15 

These  may  be  compared  to  the  results  of  a  regression 
analysis  of  the  data  in  Table  A-2 ,  based  on  Equation  41,  with 

r  *  .44  and  k  ■  5.56  x  10  4: 

0  «  0.770 
b  -  10.70 

These  results  are  plotted  in  Figure  10.  The  calculation 
was  also  carried  out  for  the  entire  fragment  population  (Tables 
A-2  through  A-5)  with  very  little  change:  P  -  .770, 

b  -  10.69  (r  -  .42,  k  -4.69  x  lO-4). 

D.  Method  IV 


This  is  the  easiest  method  to  apply,  and  potentially  the 
most  inaccurate.  It  assumes  the  relationship  between  fragment 
area  and  mass  to  be  the  same  as  that  for  spheres,  and  measures 
no  fragment  areas.  To  summarize  the  three  simplified  methods, 
each  individual  fragment's  velocity  computation  is  based  on  a 
relationship  between  that  fragment's  area  and  its  mass 
determined  from: 

Method  I  and  II.  Direct  measurement  of  each  fragment's  area 
and  mass. 
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FIGURE  10.  Results  of  Method  III. 


Method  III.  Measurement  of  fragment  area  for  a  sample  of  the 
entire  population;  but  measurements  of  all  fragment  masses. 


Method  IV.  Measurement  of  no  fragment  areas,  but  measurement 
of  all  fragment  masses;  relationship  is  that  of  spheres. 


Equation  38  forms  the  basis  of  the  approximation  used  to 
avoid  the  need  of  measuring  each  fragment's  presented  area. 

A  =  kMr 

k  =  A/Mr 
For  spheres: 


The  dependence  on  the  radius  will  be  removed  if  r  =  2/3«  In 
words,  the  proportionality  between  presented  area  and  mass  of  a 
sphere  will  be  Independent  of  the  size  of  the  sphere.  This  is 
assumed  to  hold  for  all  classes  of  homologous  fragments.  (See 
also  reference  6  for  a  detailed  discussion  of  shape  factors). 


For  steel  spheres: 
k  -  <7833,-^  (|)  2/30'/3  (r)° 
k  =  30.65  X  10_J| 

Again  substitute  Equation  38  into  Equation  25,  set  r 
rearrange : 


(42) 
2/3  and 


(43) 


(44  ) 


6.  Dehn,  James  T.  "Terminal  Effectiveness,  Vulnerability 
Methodology  and  Fragmentation  Warhead  Optimization.  I.  A  Survey 
from  an  Historical  Perspective."  ARBRL-TR-02234  ,  April  1980  . 
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Here  b  is  defined  in  Equation  37,  the  values  for  0  and  c  were  found  with  the 
sphere  calibration  data,  and  the  value  for  k  is  that  of  Equation  42.  Again, 
regression  analysis  using  measured  velocity  data  determines  the  experimental 
values  for  0  and  b.  Equation  44  is  the  basis  this  time.  The  result  is: 

0  -  0.763 

b  -  11.02 

The  data  and  regression  fit  are  plotted  in  Figure  11  ("actual")  along  with 
the  "predicted"  fit  (0  -  .672,  b  -  10.15).  This  is  a  remarkably  good  fit 
considering  the  arbitrary  nature  of  the  r  -  2/3  assumption. 


VII.  CONCLUSIONS 

Table  5  summarizes  the  various  methods  discussed,  as  applied  to 
Collins'  data  (ref  4).  It  shows  the  values  of  the  empirical  constants 
("calibration  values"),  which  are  computed  from  the  calibration  data,  and  it 
shows  the  equations  in  which  they  are  used. 

The  slopes  and  intercepts  of  the  "actual"  curves  in  figures  3,  8,  10, 
and  11  are  listed  as  the  "true  values"  in  Table  5.  These  values  were  found 
by  substituting  the  measured  data  (fragment  area,  mass,  striking  velocity 
and  penetration)  into  the  relationships  shown  on  the  axes ,  and  performing  a 
linear  regression  of  the  results. 

Two  techniques  are  used  to  compare  the  accuracy  of  the  methods .  The 
first  is  based  on  visualization.  In  a  graphical  approach,  the  actual  and 
predicted  curves  in  one  of  the  figures  (3,  8,  10,  or  11)  are  compared 
visually  to  assess  the  accuracy  of  the  method.  Similarly,  an  evaluation  can 
be  made  using  the  calibration  values  and  the  true  values  found  in  Table  5. 
For  example,  Method  II  is  more  accurate  than  Methods  III  or  IV  because  the 
empirical  (calibration)  values  for  f)  and  b  are  closest  to  the  true  values 
(those  that  result  from  using  measured  velocities  in  the  regression  analysis 
instead  of  velocities  computed  with  the  equations  listed  in  Table  5) . 

For  comparison  purposes,  Table  5  includes  information  about  the  method 
used  by  Collins  (and  others).  The  "calibrated  values"  (empirical  constants) 
used  in  this  method  are  discussed  on  page  32  (Method  II).  The  "true  values" 
given  for  (3  and  b  are  actually  the  result  of  a  linear  regression  of  the 
velocities  computed  by  the  equation  given  in  the  table.  That  is,  they  are 
not  based  on  measured  velocities,  as  are  the  other  values  given  in  these  two 
columns.  They  are  included  here  for  comparison  only;  the  main  purpose  for 
including  this  last  line  in  Table  5  is  for  the  accuracy  comparison,  which  is 
discussed  below. 

The  second  technique  for  comparing  the  accuracy  of  the  various  methods 
is  based  on  the  linear  regressions  as  well,  but  is  more  easily  related  to  an 
experiment.  This  approach  compares  the  range  of  true  velocities  for  which 
the  given  model  predicts  velocities  within  10%  of  the  true  velocity.  It  is 
based  on  the  regression  of  each  model's  predictions  as  compared  to  that  of 
the  actual  data.  These  values  are  presented  in  the  last  two  columns  of 
Table  5.  One  immediately  sees  that  the  easier  to  apply  models  are  less 
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of  Method  IV. 


accurate.  The  model  studied  in  Method  I  results  in  the  widest  accurate 
range.  As  expected  (see  discussion  of  the  technique  above),  it  is  the  best 
for  the  higher  velocities.  This  is  due  to  the  influence  of  the  exponent  of 
v  in  the  resistance  force  relationship  (Equation  1),  and  the  fact  that,  at 
high  velocities,  the  exponent  is  most  accurately  described  in  this  method. 

Method  II  sacrifices  accuracy  at  high  velocity  by  the  simplifying 
assumption  of  a  constant  exponent  of  velocity  in  the  resistance  force 
expression  (Equation  21).  But  it  is  a  more  accurate  description  at  the  low 
velocities,  and  it  is  an  acceptable  alternative  in  the  current  work. 

As  discussed  above,  Method  III  can  be  nearly  as  accurate  as  Method  II 
if  properly  applied,  and  the  accuracy  figures  support  this  contention. 

Method  IV  has  a  poor  accuracy  range,  and  at  best  should  only  be  used 
for  the  very  low  velocity  fragments.  However,  it  does  provide  some  insight 
into  the  type  of  fragments  that  came  off  of  this  weapon  at  low  velocities, 
in  that  the  shape  factor  used  was  reasonable  for  these  fragments. 

The  equation  used  by  Collins  in  reference  4  (discussed  in  section  IV 
above),  has  an  accuracy  of  slightly  less  than  Method  III.  Method  II,  with 
fewer  approximations  than  Method  III,  shows  a  greater  advantage  in  accuracy. 
The  fact  that  Method  II  considers  all  exponents  to  be  equal  in  magnitude 
lends  credence  to  the  discussion  on  that  subject  found  in  section  IV. 

The  conclusions  drawn  above  apply  to  the  use  of  the  equations  with  one 
particular  set  of  data,  and  do  not  necessarily  hold  for  other  data  sets. 

Before  closing,  two  additional  approaches  are  discussed  very  briefly. 
First,  the  possibility  of  using  a  Taylor  Series  expansion  of  f  in  Equation  1 
seems  to  have  merit,  but  has  a  major  problem:  the  derivative  of  f  does  not 
simplify  the  integral  in  equation  4,  since  (eq.  4.2.53,  ref  7): 


d_  v  ,,.i  ,  v 

~r~  v  -  (1  +  In  v)  v 


Second,  the  following  empirical  formulation  for  a  velocity  dependent 
resistive  force  was  suggested  by  J.  Dehn  as  a  possible  alternative. 


_ fc-Cy/Yo.) _ 

ln[l+(v/v0) ] 


L'Hospital's  rule  shows  that  this  force  has  a  limiting  value  of  k  as  v  -+  0, 
and  a  limiting  value  of  infinity  as  v  ■*  «,  These  limits  are  the  same  as 
those  for  the  Poncelet  equation  (f  -  a  +  bv  +  cv2  +  ...).  Substitution  into 
Newton's  law  yields  the  following  (the  same  as  Equation  25,  with  0  -  h) . 
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APPENDIX 


The  following  data  is  from  experiments 
reported  in  reference  4. 


conducted  by  Collins 


Table  A-1.  Calibration  of  Impact  Medium  using  57mg  Steel 

Spheres . 


Penetrati on 


Measured 


Penetration  j 


Measured 
Velocity  ! 


5 .  3 


6.09 

5.83 
5.30 
5.96 

5.96 
5.30 
5.56 
5.25 
5.43 
5.30 
5.30 
5.30 

3.97 
4.19 
4.19 
4.19 
3.88 

3.84 
3.84 
3.84 
3.97 
4.11 


0.14 
6 
6 

0.16 

0.15 

0.13 

0.15 

0.15 

0.13 

0.14 

0.13 

0.14 

0.13 

0.13 

0.13 

0.10 

0.11 

0.11 

0.11 

0.10 

0.10 

0.10 

0.10 

0.10 

0.01 

0.06 


3.0 

5 

0  .0 

8 

1.72 

0.04 

5036 


5256 

5211 

50  1  8 
5140 
4395 
4432 
4510 
4493 
4  367 
444  1 
4468 
3257 
3330 
3357 
3209 
3304 
3180 
3156 
3156 
324  3 
3263 
2025 
2583 


1602  i 
1  588  I 
1479 
1529 
1567 
1340 
1  351 
1375 
1369 
1331 
1354 
1362 
993 

1015  I 
1023  i 
978  i 
1007 
969 
962 
962 
988 
955 
617 
787 


3.05 
2.78 
1.85 
2.52 
1  .99 
2.25 
2.25 


2.78 
3.05 
2.78 
2.52 
1 .85 
2.52 
2.52 
1.33 
0.66 
0.66 
0.93 
0.93 
0.66 
0.93 
0.40 
0.40 
0.27 
0.27 
0.13 
0.13 


0.08 
0.07 
C  .05 
0.06 
0.05 
0.06 
0.06 
0.07 
0.08 
0.07 
0.06 
0.05 
0.06 
0.06 
0.03 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.01 
0.01 
0.0  1 
0.01 
0.003 
0.003 


2469 

753 

2139 

652 

1678 

511 

2204 

672 

1860 

567 

1853 

565 

2017 

615 

2174 

663 

2275 

693 

2266 

691 

2007 

612 

1762 

537 

2039 

621 

2139 

652 

1275 

389 

970 

296 

926 

282 

1208 

368 

961 

293 

881 

269 

1127 

343 

1093 

333 

709 

216 

661 

201 

478 

146 

■  1  ■ "  - .  J  . .  *— — - 

:  Table  A-2 .  Small  Fragment  Data:  36mg  <  m  <  80mg. 


measured 


mass,  M 

area , 

A 

!  penetration,  P 

i  velocity 

*  vo 

-5 

_  CL  O 

-6  2 

f 

1 

10  kg 

*4 

o 

c* 

10  m 

!  in 

* 

m 

1  ft/s 

m/s 

6.10 

8.81 

8.18 

1  0.93 

.024 

| 

1  1333 

406 

6.60 

7.74 

7.19 

\  0.40 

.010 

1  677 

206 

7  .00 

8.00 

7.43 

i  0.27 

.007 

!  485 

14 

6.60 

7.98 

7.4  1 

1  0.13 

.003 

!  299 

91 

5.80 

7.72 

7.17  . 

}  0.13 

.003 

!  297 

91 

3.60 

7.22 

6.71 

I  0.13 

.003 

!  219 

67 

7.50 

9.63 

8.95 

1  4.77 

.121 

i  527  2 

1607 

8.00 

10.96 

10.18 

!  4.37 

.111 

!  4645 

14  15 

7  .00 

9.15 

8.50 

i  3.31 

.084 

!  4816 

1468 

7.00 

8.99 

8.35 

!  4.50 

.114 

1  5071 

154  6 

7.00 

9.11 

8.46 

i  4.50 

.114 

!  4928 

1502 

7.00 

8.86 

8.23 

!  5.03 

.128 

!  5220 

1591 

7.00 

9.11 

8.46 

!  3.58 

.091 

I  4465 

1361 

7.10 

8.89 

8.26 

i  4.50 

.114 

!  5095 

1553 

7.10 

9.27 

8.61 

!  2.39 

.061 

!  4665 

1422 

7.10 

9.23 

8.57 

!  3.44 

.087 

!  4983 

1824 

7.10 

8.93 

8.30 

i  3.97 

.101 

!  5129 

1563 

7.10 

9.13 

8.48 

i  3.44 

.087 

!  5289 

1612 

7.10 

9.32 

8.66 

I  2.92 

.074 

i  5022 

1531 

7.10 

9.52 

8.84 

!  3.05 

.077 

!  5128 

1563 

7.20 

9.26 

8.60 

i  2.78 

.071 

!  3958 

1206 

7.30 

9.13 

8.48 

!  3.84 

.098 

!  4056 

1236 

7.40 

9.63 

8.95 

1  2.39 

.061 

!  3942 

1201 

7.40 

9.79 

9.09 

!  2.52 

.064 

i  3951 

1204 

7.40 

9.58 

8.90 

i  2.39 

.061 

i  4  17  2 

1272 

7.50 

9.06 

8.42 

!  2.52 

.064 

i  4  224 

1272 

7.50 

9.13 

8.48 

!  2.39 

.061 

!  3991 

1216 

7.60 

9.38 

8.71 

!  2.39 

.061 

!  3893 

1187 

7.60 

9.56 

8.88 

!  3.44 

.087 

!  4006 

1221 

7.60 

9.52 

8.84 

i  1.99 

.051 

i  2774 

845 

7.80 

9.63 

8.95 

3.05 

.077 

i  3165 

965 

6.90 

9.18 

8.53 

i  0.80 

.020 

i  1756 

535 

7.00 

9.02 

8.38 

i  1.99 

.051 

i  2120 

646 

7.00 

9.28 

8.62 

S  1 .46 

.037 

1  1920 

585 

7.00 

9.18 

8.53 

!  1.85 

.047 

1  1987 

606 

6.90 

9.06 

8.42 

!  1.99 

.051 

!  2232 

680 

6.90 

9.46 

8.79 

1  1.99 

.051 

i  1925 

587 

5.80 

8.90 

8.27 

}  1.46 

.037 

!  1861 

567 

6.50 

8.81 

8.18 

1  1.19 

.030 

!  1311 

400 

6.60 

8.93 

8.30 

S  1.19 

.030 

!  1292 

394 

6.70 

9.18 

8.53 

!  0.13 

.003 

i  129  2 

394 

6.80 

8.93 

8.30 

\  0.27 

.007 

1  778 

237 

6.80 

8.82 

8.19 

l  0.40... 

.010 

1  893 

272 

40 


Table  A-2 


Small  Fragment  Data:  36mg  <  m  <  80mg  (continued) 


i  6.90 

!  9.23 

8.57 

1.19 

.030 

!  1690 

515  i 

:  6.90 

!  9.05 

8.41 

0.93 

.024 

i  1577 

481  : 

!  6.90 

!  9.59 

8.91 

1  .46 

.037 

i  2093 

638  I 

!  6.70 

1  8.63 

8.02 

0.13 

.003 

!  466 

142  : 

i  6.70 

I  8.83 

8.20 

0.27 

.007 

!  515 

157  S 

1  6.70 

!  8.73 

8.11 

0.13 

.003 

!  431 

131  ! 

!  6.80 

!  8.53 

7.92 

0.13 

.003 

i  539 

164  ; 

I  _ 

_  ■ _  » _ 

Table  A-3.  Calibration  of  Impact  Medium  using  700mg  Steel 

Spheres . 


Measured  Measured 


!  Penetration 

!  Veloci ty 

Penetration 

Velocity  ! 

\  in 

m 

f  t  /s 

d/s 

i  n 

D 

!  f  t/s 

D/S  J 

700mg 

1 

i  6.76 

0.17 

2043 

623 

14.57 

0.37 

4306 

1312  { 

i  7.02 

0.18 

2024 

617 

16.03 

0.41 

4274 

1303  ! 

1  6.49 

0.16 

2020 

616 

15.77 

0.40 

4663 

1421  : 

6.23 

0.16 

1860 

567 

17.89 

0.45 

4780 

1457  ! 

i  8.35 

0.21 

2485 

757 

16.30 

0.41 

4925 

1501  i 

!  8.75 

0.22 

1976 

602 

17.09 

0.43 

4892 

1491  ! 

i  8..  6  1 

0.22 

1946 

593 

19.21 

0.49 

4950 

1509  I 

!  10.33 

0.26 

2210 

674 

17.75 

0.45 

4950 

1509  ! 

i  7.02 

0.18 

2076 

633 

19.21 

0.49 

5079 

1548  i 

!  6.76 

0 .17 

2042 

622 

19.48 

0.50 

5074 

1546  ! 

!  6.09 

0.15 

1962 

598 

14.97 

0.38 

4468 

1362  i 

!  8.88 

0.23 

2025 

616 

18.95 

0.48 

5115 

1559  ! 

!  2.12 

0.05 

997 

298 

17.62 

0.45 

5139 

1566  ! 

!  1.85 

0.05 

885 

270 

17.62 

0.45 

5018 

1529  1 

i  2.39 

0.06 

988 

301 

18.15 

0.46 

5068 

1545  ! 

i  0.80 

0.02 

621 

189 

16.30 

0.4  1 

5145 

1568  ! 

'  1.85 

0.05 

847 

258 

16.56 

0.42 

4289 

1307  I 

!  0.13 

0.003 

250 

76 

!  18.15 

0.46 

4135 

1260  ! 

i  2.78 

0.07 

1146 

349 

13.91 

0.35 

3838 

1170  ; 

S  1.59 

0.04 

624 

190 

16.83 

0.43 

3868 

1179  ! 

!  0 . 66 

0.02 

534 

163 

15.50 

0.39 

4070 

1240  I 

:  i.i9 

0.03 

569 

173 

15.63 

0.40 

4160 

1268  ! 

i  1  .72 

0.04 

809 

247 

15.11 

0.38 

4155 

1266  | 

1  2.25 

0.06 

975 

297 

14.71 

0.37 

4037 

1230  I 

i  0.66 

0.02 

571 

174 

13.91 

0.35 

4001 

1219  i 

:  2.92 

0.07 

1150 

351 

14.71 

0.37 

4004 

1220  i 

i  0.93 

0.02 

711 

217 

14.97 

0.38 

3994 

1217  ! 

!  0.93 

0.02 

647 

197 

14.44 

0.37 

3246 

989  ! 

!  2.25 

0.06 

896 

273 

14.70 

0.37 

3211 

979  ! 

!  1.33 

0.03 

764 

233 

10.87 

0.28 

2959 

902  ! 

!  0.80 

0.02 

507 

155  f 

14.57 

0.37 

3131 

954  ! 

669  «  m  <  7 OOmg 

10.87 

0.28 

3057 

932  ! 

i  12.06 

0.31 

3034 

925 

12.06 

0.31 

3106 

947  : 

1  10.33 

0.26 

3025 

922 

1 1  .79 

0.30 

2990 

911  ! 

!  9.67 

0.25 

2282 

696 

10.33 

0.26 

1961 

902  ; 

42 


Table  A-4 .  Calibration  of  Impact  Medium  using  1025mg  Steel 

Spheres . 


Measured  Measured 


.!  Penetration 

Velocity 

t 

.1 

Penetration 

I 

Velocity  i 

1  i  n 

j 

b> 

I  f  t/s 

ffl/8 

1 

1 

1 

in 

m 

~r 

ft/s 

m/s  1 

;  16.16 

0.41 

4426 

1349 

i 

» 

1 

19.48 

0.50 

» 

1 

5993 

1827  ! 

!  13.38 

0.34 

3452 

1052 

1 

1 

18.68 

0.47 

1 

5628 

1715  1 

S  15.11 

0.38 

4062 

1238 

1 

1 

20.14 

0.51 

1 

« 

6240 

1902  ' 

!  13.38 

0.34 

3528 

1075 

1 

i 

19.48 

0.50 

1 

1 

6224 

1897  1 

i  14.05 

0.36 

3553 

1083 

1 

1 

20.01 

0.51 

1 

1 

6162 

1878  } 

!  13.65 

0.35 

3473 

1059 

1 

1 

17.09 

0.43 

1 

1 

5165 

1574  ! 

!  12.32 

0.31 

3196 

974 

1 

1 

20.27 

0.52 

1 

1 

6227 

1898  | 

i  15.63 

0.40 

4075 

1242 

1 

1 

19.74 

0.50 

1 

1 

6152 

1875  : 

!  15.50 

0.39 

4033 

1229 

1 

1 

19.74 

0.50 

1 

1 

6259 

1908  ! 

!  13.12 

0.33 

337  1 

1027 

i 

18.82 

0.48 

1 

1 

5481 

1671  i 

:  n.oo 

0.28 

2855 

870 

1 

1 

16.83 

0.43 

1 

1 

4609 

1405  ! 

!  11.93 

0.30 

3063 

934 

! 

18.82 

0.48 

1 

1 

5370 

1637  i 

i  14.44 

0.37 

3990 

1216 

1 

1 

19.48 

0.50 

i 

5303 

1616  i 

I  15.50 

0.39 

4220 

1286 

1 

1 

17.09 

0.43 

« 

1 
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Table  A-5.  Calibration  of  Impact  Medium  using  1002mg  to  1059mg 

Steel  Spheres. 


Measured  Measured 
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